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Abstract 

We solve a model of self-avoiding walks which allows for a site to be visited up to two times by the 
walk on the Husimi lattice. This model is inspired in the Domb-Joyce model and was proposed to 
describe the collapse transition of polymers with one-site interactions only. We consider the version 
in which immediate self-reversals of the walk are forbidden (RF model). The phase diagram we 
obtain for the grand-canonical version of the model is similar to the one found in the solution of 
the Bethe lattice, with two distinct polymerized phases, a tricritical point and a critical endpoint. 
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I. INTRODUCTION 



Although linear polymers in solution may be studied using continuous models, much 
has been learned about these systems through models of self- and mutually avoiding walks 
placed on lattices [lj]. The excluded volume constraint makes it quite difficult to treat 
these problems, even trying to answer apparently simple questions such as the number of 
walks with a given number of steps, which is relevant in series expansion approaches to 
;he problem, is a challenging and rich field of research (for a recent work in this field see 
2I). The connection of such models with the n- vector model of magnetism in the formal 
limit n — > 3| and the application of the ideas of scaling and the renormalization group to 
polymer systems {^j have also been central in the development of this area of research. 

If the polymer chain is placed in a poor solvent, as the temperature is decreased eventually 
the chain changes from an extended configuration (in which the entropy is favored) to a 
collapsed configuration (with less contact between the polymer and the solvent and thus a 
smaller energy). The temperature where this transition happens was called 9 temperature 
Although this transition may be modeled using lattice models where the solvent is 
included explicitly these models in a certain limit lead to a simpler model which has 
become the standard model for this phenomenon, in which, besides the repulsive excluded 
volume interactions, attractive interactions between monomers on first neighbor sites but not 
consecutive along the chain are introduced (ISAW model). The configurations in this model 
are self-avoiding walks whose steps link monomers located on first neighbor lattice sites. 
The collapse transition in models which are grand-canonical with respect to the number 
of monomers in the system usually appears as a tricritical point in the phase diagram. In 
the fugacity vs. temperature phase diagram, a non-polymerized phase is present in the 
region of low monomer fugacity, and a polymerized phase is stable for higher fugacities. 
The transition between these phases is of first order at low temperatures and becomes 
continuous at higher temperatures. These regimes are separated by a tricritical point, where 
the collapse transition occurs. The ISAW model has been extensively studied on the square 
lattice js], and the exact tricritical exponents of the diluted polymer model were found 
Under certain conditions, an even richer phase diagram is found, with the presence of a 
dense polymerized phase at finite monomer fugacity, where the density of empty lattice sites 
vanishes. This additional phase was found in solutions of the model on q = 4 Husimi lattices 
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with interactions between bonds |8|, |9[ of the polymer located on opposite sides of elementary 
squares of the lattice, as well as in cluster approximations of similar models on the square 
lattice [l(3|. Transfer matrix calculations on strips of finite widths support the existence of 



this phase on the square lattice as well 



111 ] . The dense phase is absent if the interactions 



are only between monomers, but even in this case there are indications that this phase will 
be present if the polymer chains are sufficiently stiff 12]. 

Usually, the introduction of interactions in the polymer models is a source of difficulty 
both in approximate solutions and in transfer-matrix approaches. Thus, a model introduced 
recently for studying the collapse transition in polymers where only one-site interactions are 
present is quite interesting [l3j. This model allows for multiple occupancy of a site by up to 
K monomers, assuming that the attractive interactions are restricted to monomers which 
occupy the same site of the lattice. One way to justify the model would be to discretize 
the original system of a polymer in a solution using a regular lattice, and choosing the size 
of the elementary cell of this lattice to be big enough to accommodate up to K monomers 
of the polymer. Also, the length of the bonds has to be larger than the size of the cells, 
so that two monomers in the same cell will never be connected by a bond. Two versions 
of the model were studied by extensive numerical simulations in [l3j on the square and on 
the cubic lattice. In the RA (immediate reversals allowed) model, there are no restrictions 
on the walks on the lattice, while in the RF model (immediate reversals forbidden) only a 
subset of the possible walks is considered: those in which the walk does not return to the 
original site immediately after reaching a new site. In the simulations done for the RF model 
on the cubic lattice for K = 3, a transition between extended and collapsed polymerized 
phases is found. 

Recently, the RA and RF models were solved on a Bethe lattice for K — 2, in order to 
compare their thermodynamic behaviors with the much studied ISAW model regarding the 
usual collapse transition [ijj] . The solution of the RF model on the Bethe lattice produces a 
phase diagram in which the polymerization transition remains continuous for small non-zero 
values of u>2, becoming of first order at higher values of this statistical weight, a behavior 
which is also found in the ISAW model. Besides the regular polymerized phase, a second 
phase is stable for sufficiently high values of u 2 and low values of uu% where only empty and 
double-occupied sites are present. 

Here the RF model is solved on a q = 2(a + 1) Husimi lattice, which is the central region 
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of a Cayley tree built with squares. There are a + 1 squares incident on each lattice site. 
The thermodynamic behavior of models on such a lattice is expected to be closer to the one 
obtained on a regular lattice with the same coordination number, and we were motivated for 
this calculation mainly for two reasons. Although the RF model on the Bethe lattice displays 
a tricritical point which may be associated to the usual collapse transition of polymers, if we 
parametrize the model such that u)\ = z, the activity of a monomer, and u>2 = z 2 u, where 
u is the Boltzmann factor associated to a pair of interacting monomers, the tricritical point 
is found at to < 1, a value that corresponds to a repulsive interaction between monomers at 
the same site. We are interested in finding out if a calculation which should lead to results 
closer to the ones on regular lattices might shift the tricritical point to the region in the 
parameter space which corresponds to attractive interactions. Also, it is of interest to find 
out if the second polymerized phase is still present in the phase diagram of the model on 
the Husimi lattice. We found that actually with a convenient parametrization the tricritical 
point is located in the physically expected region and that the second polymerized phase 
is still present in the Husimi lattice solution, although it occupies a smaller region of the 
parameter space than the one found for the Bethe lattice. 

In section [III we define the model in more detail and present its solution on the Husimi 
lattice. Final discussions and the conclusion may be found in section IIHI 

II. DEFINITION OF THE MODEL AND SOLUTION ON THE HUSIMI LATTICE 

We consider a Husimi tree, a Cayley tree built with polygons, which in our case will be 
squares. Sometimes this tree is also called a cactus. As also happens for the Cayley tree, 
in the thermodynamic limit the fraction of sites which are on the surface of the tree does 
not vanish, and this accounts for the fact that the solution of models on such trees usually 
shows a behavior which does not resemble the one found on regular lattices. If, however, 
the behavior of models in the central region of the trees is considered, for many models the 
exact solution corresponds to the Bethe approximation of the same model on a regular lattice 
with the same coordination number, and this is the reason why this is called a Bethe lattice 
solution [151 ] . The Husimi lattice corresponds to the central region of a Cayley tree built 
with polygons (squares in our case), and since closed paths are present (although restricted 
to single elementary squares), it is expected that the solution of models on this tree will be 
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closer to the one found on regular lattices, and this is confirmed in many cases. 

The allowed configurations of the K = 2 RF model are walks that may visit a site one 
or two times, with their initial and final monomers placed only on the surface of the tree. 
Also, as stated above, immediate reversals of the walk are forbidden. In figure [T] a possible 
configuration is shown on a tree with two generations of squares. The statistical weight of 
a configuration will be where N\ and N% are the numbers of sites with one and 

two monomers, respectively. For simplicity, the surface of the tree was chosen to be defined 
by sites connected to a single site of the first generation of squares. To solve the model 
on the Husimi tree, we consider rooted subtrees, and define partial partition functions for 
these trees, for fixed configurations of the bonds incident on the root site. The operation of 
attaching three sets of a n-generations subtrees to a new root square will result in a n + 1- 
generations subtree, leading to recursion relations for the partial partition functions. Notice 
that we decided to consider the monomers placed on the same site to be indistinguishable, 
in opposition to what was adopted in the solution on the Bethe lattice [jjj] , where they were 
supposed to be distinguishable. This new convention was motivated mainly by two aspects: 
as already observed in the discussion of the Bethe lattice solution, the tricritical point of 
the K = 2 RF model is located in the region in the parameter space which corresponds 
to repulsive interactions between monomers. If two monomers placed on the same site are 
considered to be indistinguishable, the weight uj 2 will be multiplyed by a factor of 2, and 
this will shift the tricritical point towards the region of attractive monomers. Also, in 

n 

the original simulations [13J the configurations of a walk on the lattice are labeled by the 
sequence of sites visited by the walk, and this convention corresponds to indistinguishable 
monomers. Although we do no detailed comparisons of our results with the simulations, 
since they were done for K = 3, we decided to adopt the same convention. In figure [2]the 11 
configurations of the bonds incident on the root site of subtrees are depicted. For brevity, 
only one configuration of pairs related by reflection symmetry is shown. For example, there 
are a total of 4 configurations with the same partial partition function in group 7 of the 
figure. We notice that the connection of the incoming bonds to the monomers placed on the 
root site is not fixed. Since we assume the monomers to be indistinguishable, to write down 
the recursion relations for the partial partition functions we must separate incoming double 
bonds in two groups: if both bonds visited the same sites since the boundary of the tree they 
are labeled as i (as in configuration 5), otherwise, they are labeled as d (as in configuration 
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FIG. 1: Example of a configuration for the K = 2 RF model on a Husimi tree with square 
ramification a = 1 and 2 generations. Polymer bonds (steps of the walks) are represented by thick 
lines, while the lattice bonds are thin lines. The weight of this configuration is u>\ 2 u>2- 



4). This information is essential to allow us to correctly determine the multiplicity of the 
contributions to the recursion relations below. 
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FIG. 2: Root configurations of the RF model with K = 2 on the Husimi lattice. 



We may now proceed obtaining the recursion relations for the partial partition functions. 
We notice that certain partial partition functions appear in the recursion relations only 
in linear combinations. Thus, g% and g% always appear as g3 + g&, and the other linear 
combinations are 2g 7 + g 8 and 2g 9 + g n . We therefore may reduce the number of independent 
variables in the recursion relations by 3. The partial partition functions usually diverge in 
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the thermodynamic limit, thus we consider the 8 ratios defined below, which often remain 
finite, and in the thermodynamic limit a phase of the system may be associated to a fixed 
point of the recursion relations for the ratios. 



Ri = -, (la) 



R2 = -, (lb) 
9o 

D 93 + 96 /, x 

R 3 = , (lc) 

9o 

R4 = -, (Id) 
9o 

2g 7 + g$ 

R 5 = , (le) 

9o 

R G = ^±^1, (if) 



Ri = -, (ig) 

9o 

Rs = ^. (ih) 

9o 

It is now convenient to define first some linear combinations of partial partition functions 
which appear repeatedly in the recursion relations, whose contributions are shown graphi- 



7 



cally in figure [3] for the particular ramification of squares a = 1. They are: 
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FIG. 3: Contributions to the vertex functions. The monomers are represented by dots. 



Now we may proceed considering the operation of attaching three sets of a subtrees to a 
new root square, summing all possible contributions for a fixed configuration of the bonds 
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FIG. 4: Contributions to the recursion relation for the partial partition function g$.. 

which are incident on the new root site. In figure H] the graphical representation of the 
contributions to g' are shown to illustrate this process, in the order they appear in the 
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equation below. The recursion relations are: 

g' = A 3 + 2AB 2 + B 2 (2C + H + D + 2E) + (2 A + G)(H 2 + 2HC + 2C 2 + ACE + 2EH) 



+2BF(2C + H + 2E); (3a) 
g[ = 2A 2 B + 2B Z + 2AB(2C + H + D + 2E) + 2B(2C + H + D + 2E) 2 

+2AF(2C + H + 2E) + 2FG(2C + H + 2E) + 2B(H 2 + 2#C + 2C 2 + ACE + 2£if) 

+2F(2C + H + D + 2E) (2C + H + 2E) + 4RF 2 ; (3b) 
# 2 = AB 2 + 25 2 (2C + if + D + 2£) + (2C + ff ) 3 + 3(2C + H) 2 (D + 2£) 
+3(2C + H)(D + 2£) 2 + 2BF(2C + H + 2E) + F 2 G + AF 2 (2C + ff) 

+2F 2 ( J D + 2£); (3c) 

# 3 = (D + 2E) 3 + 2F 2 (D + 2E) + F 2 G; (3d) 
^ = 2A 2 C + 2ABF + 2 J B 2 C + 2BF(2C + H + D + 2E) + 2ACG + 2CG 2 + 

2F 2 (2C + H + 2E) + 2BFG + 2C(H 2 + 2HC + 2C 2 + ACE + 2£ff); (3e) 

g' h = 2H{A 2 + B 2 + AG + G 2 + H 2 + 2HC + 2C 2 + ACE + 2EH); (3f) 
g' 6 = 2A 2 E + 2B 2 E + 2AEG + 2EG 2 + 

2£(ff 2 + 2HC + 2C 2 + ACE + 2£ff ) ; (3g) 
g' 7 = 2AB(C + H + E)+ 2B 2 F + 2B(C + H + E) (2C + H + D + 2E) + 2 J BG(C + ff + E) + 
2F(C + H + £) (2C + ff + 2£) + 2F[(2C + H) 2 + 2D(2C + ff ) + 4£(2C + ff)] 

+2FG(2C + ff) + 2F 3 ; (3h) 

g' 8 = 2F(D + 2E) 2 + 2FG(D + 2£) + 2FG 2 + 2F 3 ; (3i) 

g ' g = AC(C + 2E) + 2BF[C + E) + 2CG[C + IE) + F 2 (2C + ff); (3j) 

g' w = AH(2C + H + 2E)+ 2BFH + 2GH(2C + H + 2E); (3k) 

g ' u = F 2 (D + 2E + 2G). (31) 



In these recursion relations the need to distinguish between the cases i and d of double 
bonds in the definition of the partial partition functions is apparent, for example, in the 
contributions to g' 2 : the terms proportional to C (case d) have an additional factor 2 as 
compared to the terms proportional to H (case i). 

Finally, we will consider the operation of attaching a + 1 subtrees to the central site of 
the lattice. This operation is similar to the ones realized to obtain the recursion relations 
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and the result is: 
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It is now easy to obtain the probabilities of single and double occupancy of the central site. 
The results are 
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FIG. 5: (color on line) Phase diagram of the model on the Husimi lattice. Full lines are continuous 
transitions and dashed lines are first order transitions between the NP and RP (blue) phases and 
between the PP and RP (purple) phases. The tricritical point and the critical enpoint are also 
shown. 



As stated above, to study the thermodynamic behavior of the model we have to find the 
fixed point of the recursion relations, given the statistical weights uj\ and oj 2 . As expected, 
the ratios R-j and i?§ vanish for all values of the statistical weights. On the Husimi lattice a 
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double chain which starts on the boundary may end at any step, by entering into a square 
and circulating it (this corresponds to the contribution 3HD 2 in the recursion relation [3c] 
for g 2 ). At low values of the statistical weights, a non-polymerized phase, characterized by 
Pi — P2 — 0, is stable. At the fixed point associated to this phase all ratios vanish, with the 
exception of R3. From the recursion relations above we may find the following equation for 
this fixed point: 

{u x + 2au2R 3 ) 3 + [2auj 2 {\ + uj 2 + uj\) - 1}R 3 = 0. (7) 

At the stability limit of the NP phase the largest eigenvalue of the jacobian of the recursion 
relations is equal to one, and this condition allowed us to obtain the region of the parameter 
space where the NP phase is stable. For large values of the statistical weights, a regular 
polymerized (RP) phase is stable, in which the ratios R±, R 2 , ■ ■ ■ , R& are non- vanishing at 
the fixed point. Finally, for small values of ui, between the regions where the NP and RP 
phases are stable, a third phase appears, for which all ratios which correspond to an odd 
number of incoming bonds at the root site -R5) vanish. We will call this phase PP (pair 
polymerized). In figure [5] the phase diagram of the model is presented for a lattice with 
q = 4, and we notice that the transition between the NP and RP phases may be of first or 
second order, with a tricritical point located at u)\ = 0.3325510(6) and u 2 = 0.120544(4). 
The region where the PP phase is the most stable one is rather small and the transition 
between the NP and PP phases is continuous. This transition line ends at a critical endpoint, 
which is located at oj\ = 0.0695605(5) and u 2 = 0.3370740(2). A discontinuous transition 
separates the two polymerized phases. The first order lines were obtained directly from the 
recursion relations, starting the iterations with 'natural' initial conditions [9]. In the present 
calculations, we considered the surface of the tree to be formed by sites connected to a single 
site of the next generations, as shown in figure [TJ This choice leads to the following initial 



12 



values: 
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As mentioned above, the ratios R-j and R$ vanish in the termodynamic limit, but have 
nonzero values for finite lattices. If we change slightly the initial conditions, assuming 
the monomers placed on the surface of the tree to be distinguishable, these ratios will 
vanish identically and the initial value for R4 would be equal to 2u)2, with no change in the 
thermodynamic properties of the model within our numerical precision. 

A model which is very similar to the one we are studying here was investigated recently 
by R. A. Zara and M. Pretti Q to study the properties of RNA-like molecules, and actually 
the phase diagram for the Husimi lattice solution of this model is similar to the one we obtain 
here. It may also be mentioned that no dense phase, as the ones found in some versions if 
the ISAW model , is stable in finite regions of the parameter space. Actually, such a 
phase is stable in a region of the u>i = line, but this fixed point is never reached if U\ 7^ 
and therefore is of no physical relevance. 



III. FINAL DISCUSSIONS AND CONCLUSION 



Qualitatively, the phase diagram presented here is similar to the one found for the RF 
model on the Bethe lattice (figure 3 of reference 14J )• We notice that on the Bethe lattice 
the second order line between the NP and the RP phases is located at io\ = 1/3, while 
the NP-PP transition happens at u)i — 1/6 (00% = 1/3 if the monomers are considered to 
be indistinguishable). Those lines are no longer parallel to one of the axes for the Husimi 
lattice solution. The location of the TCP is not changed much in the two solutions, although 
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the value of 002 for this point on the Husimi lattice solution is about 10 % larger. The CEP, 
although localized at almost the same value of UJ2 in the two solutions, shows a much smaller 
value of uj\ on the Husimi lattice. As a consequence, the area in which the PP phase is stable 
in the parameter space is much smaller on the Husimi lattice solution, it is an open question 
if this unusual phase appears if the model is considered on regular lattices. On the Bethe 
lattice solution the PP phase was called double occupancy polymerized phase, since p\ was 
found to vanish in this phase. We notice that in the Husimi lattice solution the density of 
sites occupies by a single monomer does not vanish in this phase, as may be appreciated in 
the inset of the figure El where both densities are shown as functions of 0J2 for a fixed value 

of UJ\. 
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FIG. 6: Density of sites occupied by one (pi, dashed lines) and two (p2-full lines) monomers, as 
functions of LO2 for uj\ = 0.05. 



The comparison of the model with multiple monomers per site (MMS) with the usual 
ISAW model is not straightforward. When 002 = 0, the MMS model corresponds to the 
ISAW model without attractive interactions. However, the MMS model for K = 2 without 
attractive interactions corresponds to the line u 2 = u\, since we should associate a statistical 
weight equal to the activity z to each monomer placed on the lattice. In the ISAW model 
a subset of the walks considered in the MMS model is allowed. If we actually imagine this 
model to be an effective description of a continuous model treated in a cell approximation, we 
might think the parameter uj to be the effective interaction found integrating the position 
of the two monomers inside the cell they occupy. If this interpretation is adopted, the 
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statistical weight of a site occupied by one monomer would be oo\ = z, where z is the 
activity of a monomer, while a double occupied site would contribute with a factor ojz 2 to 
the grand canonical partition function, where u is the Boltzmann factor associated to the 
interaction of two monomers. With this parametrization, the tricritical point will be located 
at uj pa 1.09, which corresponds to attractive interactions between monomers at the same 
site. In the Bethe lattice solution for indistinguishable monomers, the tricritical point is 
located at uj — 1, which corresponds to no interaction. 

We also did the calculations of the model for distinguishable monomers, as was done 
initially for the Bethe lattice. Within our numerical approximation, this solution leads to 
a phase diagram which differs from the one presented here by a factor of 2 in the values 
of u) 2 . This may be understood by noting that the recursion relations for the model with 
distinguishable monomers are the same ones presented here except for this factor 2 in each 
term proportional to u 2 for Ri, R 2 , . . . , R$ (configurations 5 and 10 are absent in this case), 
and that the additional ratios R-j and R$ vanish in the fixed point for the model with 
indistinguishable monomers, as stated above. 
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